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INTRODUCTION 


There  are  many  facets  to  the  problem  of  turbulent  combustion.  While  many  aspects 
of  the  problem  can  now  be  handled  satisfactorily,  a  remaining  difficulty  is  the  microscale 
processes.  These  can  be  viewed  in  different  ways.  In  terms  of  the  fundamental 
conservation  equations,  the  microscale  processes  are  mixing  and  reaction.  The  essence  of 
the  problem  is  that  (in  the  appropriate  regime)  these  processes  are  strongly  coupled  —  rapid 
reaction  being  intimately  connected  with  rapid  mixing,  for  example.  In  the  alternative 
view,  combustion  occurs  on  thin  flame  sheets.  Then  the  issues  are  to  understand  the 
statistics  and  dynamics  of  these  sheets.  The  sheets  are  convected,  bent  and  stretched  by  the 
turbulence,  and  their  local  propagation  speed  is  influenced  by  the  local  straining  and 
curvature. 

In  this  study,  Direct  Numerical  Simulations  (DNS)  of  turbulence  are  used  to  study 
basic  microscale  processes.  Even  with  current  supercomputers,  the  parameter-range 
accessible  to  DNS  is  quite  small.  If  an  inhomogeneous  flow  (e.g.  a  flame  or  mixing  layer) 
is  simulated  by  solving  the  basic  conservation  equations,  then  the  practically-interesting 
parameter  range  cannot  be  reached.  In  particular  a  good  separation  between  micro  and 
macro  scales  cannot  be  achieved,  while  (at  the  same  time)  maintaining  a  large  Karlovitz 
number.  Instead  we  study  inert  homogeneous  turbulence  and  extract  information  directly 
relevant  to  turbulent  combustion. 

In  previous  work  (Yeung,  Girimaji  &  Pope  1989,  Pope,  Yeung  &  Girimaji  1989) 
we  have  studied  the  straining  and  curvature  of  material  surfaces  —  which,  iijjSippropriate 
limits,  are  statistically  identical  to  flame  sheets.  The  accomplishments  in  this  reporting 
period  —  described  below  —  are  in  four  areas.  The  first  two  provide  more  information  on 
material  surfaces  and  related  quantities;  the  third  is  an  extension  of  these  studies  to 
propagating  surfaces;  and  the  fourth,  separate,  contribution  is  concerned  with  molecular 


mixing. 


ACCOMPLISHMENTS 
Material  Element  Deformation 

In  Girimaji  &  Pope  (1989a)  the  previous  studies  of  material  surface  elements  is 
extended  in  two  ways:  first,  additional  important  statistics  are  reported;  and,  second, 
material  line  and  volume  elements  are  also  considered.  The  major  findings  are  now 
itemized. 

1.  The  growth  rates  of  line  and  surface  elements  are  significantly  less  (by  a  factor  of 
about  3)  than  previous  estimates  (Monin  &  Yaglom  1975). 

2.  The  smaller  growth  rates  are  due  to  the  poor  alignment  between  material  elements  and 
the  turbulent  straining  motions.  The  lack  of  alignment  is  caused  both  by  the  action  of 
vorticity  to  rotate  the  elements,  and  because  the  principal  axes  of  the  strain-rate  rotate 
quite  rapidly.  In  other  words,  contrar_,  to  conventional  wisdom,  we  find  that  strain  is 
fleeting  rather  than  persistent. 

3.  An  initially  spherical  infinitesimal  volume  element  is  deformed  by  the  turbulent 
straining  into  an  ellipsoid.  The  most  probable  shape  of  this  ellipsoid  is  like  a  squashed 
cigar  —  one  axis  is  extended,  another  equally  compressed,  while  the  third  remains 
(approximately)  unchanged. 

Stochastic  Model  for  Velocity  Gradients 

Any  quantity  related  to  material  element  deformation  —  the  straining  on  a  material 
surface,  for  example  —  can  be  calculated  from  the  velocity-gradient  time  series  following  a 
fluid  particle.  Thus  the  study  of  material  element  deformation  described  in  the  previous 
subsection  was  performed  using  the  velocity-gradient  time  series  obtained  from  DNS. 

We  have  devised  a  computationally-simple  stochastic  model  for  the  velocity- 
gradient  time  series  that  reproduces  the  major  statistics  of  the  DNS  time  series.  This 
stochastic  model  can  be  used  in  turbulent  combustion  models  (e.g.  that  of  Pope  &  Chen 
1988)  to  calculate  the  strain  rate  on  the  flame  sheet  etc. 
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The  stochastic  model  (fully  described  and  tested  by  Girimaji  &  Pope  1989b)  is  a 
tensor-valued  diffusion  process.  Thus  it  is  a  Markov  process  with  continuous  sample 
paths.  The  drift  and  diffusion  coefficients  in  the  model  have  been  determined  —  partly 
analytically,  partly  empirically  —  by  requiring  that  the  most  important  statistics  of  the 
process  match  those  obtained  from  DNS.  Typically  the  model  is  accurate  to  within  15%. 
Propagating  surfaces 

The  material  surfaces  described  so  far  represent  flame  sheets  only  in  special  limits. 
In  the  context  of  premixed  flames,  the  limit  is  ul/ut|— >0,  where  ul  is  the  laminar  flame 
speed  and  u^  is  the  Kolmogorov  vdov-dy  scale. 

We  have  nearly  completed  a  study  of  propagating  surfaces,  with  propagation 
speeds  ul  of  u\J\ir\  =  0,  j,  1,  4,  16.  The  study  of  propagating  surfaces  (ul>0)  is  much 

more  difficult  than  that  of  material  surfaces  (ul=0),  because  of  cusp  formation.  That  is.  the 
curvature  of  a  surface  element  can  become  infinite  in  finite  time. 

The  first  difficulty  to  be  faced  is  to  devise  an  accurate  and  stable  numerical 
algorithm  that  can  calculate  the  surface  properties  up  to  the  point  of  cusp  formation.  This 
requires  special  consideration  of  the  singularity  (i.e.  the  cusps).  Such  an  algorithm  has 
been  devised.  It  is  second-order  accurate,  and  reduces  to  an  exact  analytic  solution  as  the 
singularity  is  approached. 

Because  of  the  formation  of  cusps,  the  time  series  of  surface  properties  are  finite  in 
duration,  and  they  are  inherently  non-stationary.  Thus  the  second  difficulty  is  the  analysis 
of  these  time  series  —  as  compared  to  the  -  M  ,iary  time  series  of  arbitrarily  long  duration 
for  material  surfaces. 

The  preliminary  results  confirm  expectations.  For  example,  the  Lagrangian  time 
scales  and  the  straining  on  the  surface  decrease  with  increasing  UL/urp  The  precise 
variation  of  these  statistics  —  yet  to  be  quantified  —  is  of  fundamental  importance  to 
theoretical  questions  of  premixed  flame  propagation.  Also  as  expected,  the  mean  time  to 
cusp  formation  decreases  with  increasing  ul/ut|- 


The  major  effort  for  the  first  6  months  of  1990  will  be  to  fnalize  and  report  the 
results  of  this  study  of  propagating  surfaces. 

Kraichnan's  PDF  Model 

For  decades  to  come,  direct  numerical  simulation  will  not  be  used  as  an  engineering 
design  tool  for  turbulent  combustion  devices,  because  the  computation  demands  are  wav 
beyond  reach.  Instead,  our  philosophy  is  to  use  DNS  to  guide  and  test  the  development  of 
statistical  models  that  can  be  used  as  design  tools.  For  turbulent  combustion  the  most 
favorable  statistical  approach  is  the  pdf  method  (see,  e.g.  Pope  1985). 

In  the  pdf  approach,  while  complex  reactions  can  be  treated  without  approximation, 
mixing  —  i.e.  molecular  diffusion  processes  —  have  to  be  modelled.  For  the  last  1 5  years 
the  development  of  a  complete'y  satisfactory  mixing  model  has  proved  an  allusive 
objective.  Recently  Kraichnan  (1989)  has  proposed  a  completely  new  modelling 
methodology  which  appears  to  be  a  breakthrough.  As  described  in  the  following 
paragraphs,  we  have  obtained  an  analytic  solution  to  Kraichnan's  model  equations  which 
are  in  spectacular  agreement  with  our  previous  DNS  results.  First  these  DNS  results  are 
reviewed. 

Eswaran  &  Pope  (1988)  (work  supported  by  AFOSR  85-0083)  performed  direct 
numerical  simulations  of  the  decay  of  an  inert  passive  scalar  0(x,t)  in  stationary, 
homogeneous,  isotropic  turbulence.  The  initial  condition  corresponds  to  blobs  of  scalar 
with  0=1  and  an  equal  quantity  of  blobs  with  0--1.  Thus  the  initial  pdf  of  0.  f(\j/), 
corresponds,  approximately,  to  the  double-delta-function  distribution. 

f(v|/)  =  ^[5(\|/+l)  +  5(\|/-1)1  .  (1) 

The  results  show  that  the  decay  of  the  scalar  variance  depends  on  the  size  of  the 
initial  blobs  (relative  to  the  turbulence  scales)  —  consistent  with  experimental  observations. 
However,  it  is  found  that  the  shapes  adopted  by  the  pdfs  are  independent  of  blob  size. 
Figure  1  shows  the  pdf  s  obtained  from  the  DNS  for  three  different  values  of  the  initial 
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blob  size.  It  may  be  seen  that  the  shapes  are  essentially  the  same  in  each  of  the  three 
figures  (Figs.  la,b,c). 

Kraichnan’s  model  is  conceptually  simple,  though  technically  difficult.  The  idea  is 
to  map  the  non-Gaussian  scalar  field  0(x,t)  to  a  Gaussian  field  <J>0(xd)  by  a  mapping  Y(*,t): 
i.e.  0o(x,t)  =  Y(4>[x,t],t).  Gaussian  closure  approximation  are  then  applied  to  the  Gaussian 
field  0O. 

For  the  case  of  an  inert  passive  scalar  with  the  initial  condition  given  by  Eq.  ( 1 )  we 
have  obtained  the  following  analytic  solution  to  Kraichnan's  model: 

f(\]/;t) = \  i  exp  { -  4  Vo  (i  -S'2) )  •  r-  - 

where 

=  el*-l  , 

Vo  =  V2  Z  erf-](\j/) , 

and  t*  is  normalized  time.  Figure  2b  shows  the  pdf  s  according  to  Eq.  (2)  compared  (Fig. 
2a)  to  those  from  DNS.  It  may  be  seen  that  the  agreement  is  remarkable.  In  spite  of  the 
considerable  efforts  of  several  workers  over  the  past  15  years,  no  other  models  have  come 
close  to  Kraichnan's  in  qualitative  or  quantitative  performance. 
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:IG.  la  :The  scalar  pdf  from  simulation  F2a  (Re.,  ~50,  kykt,  =  1 )  at  dif- 
irent  times.  Lines:  —  (tu/l  ~  0.22;  tb'/i „  =0.99),  ( 1.49;  0.73),  -  +  - 

2.11;  0.55), -o-  (2.78;  0.40), -x-  (3.47;  0.27). 


FIG.  lb  :The  scalar  pdf  from  simulation  F2c  ( Re,,  ~  50,  kjk^  =  4)  at  dil 
ferent  times.  Lines:  —  ( tu/l  =  0.22;  <t>'/<t>'0  =  0.92),  —  (0.42;  0.80),  -  + 
(0.83;  0.54),  -o-  ( 1.28;  0.35),  -x-  ( 1.49;  0.28). 


Fig.  1:  Scalar  pdf's  from 
direct  numerical  simulations 
of  Eswaran  &  Pope  (1988). 
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FIG.lc  :  The  scalar  pdf  from  simulation  F2e  (Re,,  ~50,  ks/k„  =  8)  at  dif¬ 
ferent  times.  Lines:  —  ( tu/l  —  0.07;  <f>'/<t>'n  =  0.94),  —  (0.22;  0.76),  -  -I-  - 
(0.42;  0.54),  -o-  (0.62;  0.38),  -x-  (0.83;  0.27). 
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